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We obtain the hadronic mass spectrum in the 'bag of bags' statistical bootstrap 
model (BBSBM), implementing the colorless state condition, aside of baryon and 
strangeness conservation, using group projection method. We study the partition 
function, investigate the properties of dense hadronic matter, and determine the 
conditions under which the system undergoes a phase transition to a deconfined 
quark-gluon plasma. We show that a phase transition cannot occur in the N = 1 
(Abelian) limit of our model, and is first order for QCD-like case N = 3. 



Dedicated to the memory of Peter A. Carruthers 
1 Introduction 

Considerable theoretical and experimental effort is being devoted to the study 
of the phase transition between a hadron gas and a quark-gluon plasma (QGP). 
This phase transition is believed to have occurred in the early Universe, and is 
searched for in high energy nuclear collisions in laboratory. The fundamental 
theoretical treatment of the hadronic phase transition relies on lattice QCD ap- 
proximation. This approach exhibits a phase transition and gives the equation 
of states of both phases. The case of nonzero baryon density (finite chemical 
potential) comprises not yet solved technical difficultiestl 

In a phcnomcnological study of this case one starts from the low tempera- 
ture phase, i.e., the hadron phase, and implement the high temperature phase 
in a more or less ad hoc way. This was the case of the statistical bootstrap 
model □ (SBM) when a quark component is added to ft beyond the critical 
curve limiting the existence domain of the hadron phase u. Here we show that 
it is possible to extend SBM in a natural way to encapsulate the high temper- 
ature deconfined phase. We next recall the SBM and in particular illustrate 
derivation of the hadronic mass spectrum. We introduce the modifications 
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leading to a SBM-type model which describe both hadron phase and QGP 
phase in section |^, and we study the phases and properties of the system at 
finite baryochemical potential in section || 



2 Statistical bootstrap 

2.1 The hadronic mass spectrum 

The SBM approach can be characterized in three statements: 1) hadrons are 
made of hadrons; 2) compound hadrons and 'elementary' hadrons have to be 
treated on the same footing: this is the bootstrap; 3) the main effect of strong 
interactions is hadron production. Hadron-hadron interactions in thermally 
and chemically equilibrated system can be accounted for (nearly) completely 
by considering the abundances of particles present, i.e., a strongly interact- 
ing hadron gas can be described as a mixture of infinitely many ideal gases, 
weighted by the hadronic mass spectrum. The dynamical problem in SBM is 
the theoretical determination of the hadronic mass spectrum, i.e., the number 
of hadron states per unit interval mass. 

For pointlike .particles in a box of volume V the bootstrap equation for 
this density reads □ 



In this equation, p out represents the total density of states in the box, while 
Pin(rrii) counts the number of mass states of particle i. Note that the density 
of levels is written with center of mass at rest because this is the density to be 
identified with the number of hadron states per unit interval of rest mass. 
One has to solve Eq. (fll) with the bootstrap condition: 





(1) 



pout(m) = pin(m) = p(m) ■ 



(2) 



The solution is of the form 



p(m) 



cm e 



with 



= -3. 
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2.2 The partition function 

The partition function in the Boltzmann approximation can be written 
lnZ(y,T) = g^ p(m)m 2 K 2 ^ dm, 



(4) 



where Ki is the modified Hankel function. With a mass spectrum of the 
form (|^) , one finds for T — > Tq, 
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and the energy density, 
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One sees that for a > —7/2, To is attained with infinite energy density and 
appears as a maximal limiting temperature. All energy delivered is converted 
into particle multiplicity rather than thermal motion. But for a < —7/2, T is 
reached at finite energy density and the system may undergo a phase transi- 
tion. In the original SBM a > —7/2 and thus To is a limiting temperature Bo. 
Interestingly, the hadronic partition function in SBM is singular (a square- 
root singularity) even for finite volume. This last feature can be traced back 
to the assumption of pointlike hadrons. However when the individual hadron 
volumes are taken into account through the so called 'excluded volume' correc- 
tion To is reached at finite energy density but not in the usual thermodynamic 
limitta. When baryon number is considered, there is not a single transition point 
To, but a transition curve T(/i) in the plane T — /i, where \x is the baryonic 
chemical potential, and T = T(fi = 0). Beyond this critical curve the SBM 
partition function becomes complex and extraneous assumptions are needed 
to implement the new QGP phase. This results in a two component model, 
one component describing the hadron side and another the quark-gluon side. 
However, a single closed analytic model unifying both phases can be built, 
provided that we introduce from the start the quark structure of hadrons as 
extended objects. 
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3 Bag of bags SBM 

3. 1 The bag grand canonical partition function 

Hadrons are considered now as bags of quarks and gluons occupying finite vol- 
umes. They have to be colorless objects, i.e., SU C (3) (or SU C (N)) singlets. 
This is an expression of color confinement. To be realistic, we also implement 
baryon number and strangeness conservation effects. With these constraints, 
we will determine the hadronic mass spectrum in such a bag-of-bags-SBM (BB- 
SBM). Consider a colorless bag of volume v filled with quarks, anti-quarks and 
gluons at sufficiently high temperature to neglect their interactions (asymptotic 
freedom). The density of states of this bag can be determined by taking the 
inverse Laplace transform of its grand canonical partition function (GCPF). 
But quarks and gluons transform respectively under the fundamental and the 
adjoint representations of SU C (3). In the GCPF one has to retain only the 



singlet states. Furthermore, as we have noted in Sec. 2.1, the density of levels 
is defined with the center of mass at rest. 

With the preceding constraints the GCPF, can be derived, using projection 
techniqueB, from the function 

Z Q {P,v,n q ,ns)= f d»(g) f — {Tr G U G {g)e-^+^ } (8) 

JSU C (3) Jv V I ) 

x {TrM a ^)e- pa '^'- A ^'}{Tr s U s ^)e-^ a '^-^' S '\ , 

where G stands for gluons, q for the quarks u and d and s for the strange quark. 
N q and N s are the quark number operators and fi q and fi s the corresponding 
chemical potentials. Pq , P q and R are the momenta of gluons and quarks and 
their space position inside theJpag. U(g) is a unitary representation of SU C {3) 
in the Hilbert space of states 19, g being an element of SU C (3). The physical 
grand partition function Zq((3,v, [Ag, /x s ) is obtained from Zq{/3,v, /j, q , fj, s ) by 
the so-called 'Wick rotation' /i — * — i/3/i. In spite of its apparent complexity 
Eq. (|^) corresponds simply to the product of five partitions functions (five 
traces) for non interacting five particle species (G, q, q, s and s). The integra- 
tion over the SU C (3) group with the Haar measure d^,(g) selects the SU C (3) 
singlet states and the integration over R in the volume v selects the bag states 
of zero momentum. 

Zq((3, v, [iq, /Lt s ) can be written in form 0, 

Z QG (p,v,fiq,fi s )= [ dfj,(g) [ A^e. e , (9) 

JSU C (3) Jv/P* V /P 
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where r = R/f3 and 



6 = ^3 (1 + r2)2 [ d G Q(g) + d q Q{g, ng) + d s Q(g, /*.)] , (10) 

with 

oo 1 

= £ £4X1,1(0*) > (11) 

fe=0 



for gluon contribution and 

Q(g^)^2nY / K -^^ X{m (g k ), (12) 



( _ 1)fc -i 



A: 4 

k=0 

for quark contribution. Xi,oG?) an< i Xi,i(<?) = Ixi,o(s , )| 2 ~~ 1 are respectively 
the characters of the fundamental and the adjoint representation of SU C (3). 
As a class function the characters depend on the eigenvalues of e %di of g, (i=l, 

2, 3): xMg) = ES e * with = °- In E q- ©» <fc = 2, = 4 

and <i s = 2 are the number of degrees of freedom of gluons and quarks, not 
counting the color multiplicity, already taken into account in the characters 
since 0) = 8 and Xi,o(0, 0) = 3. 

For high temperature, (/? — > 0), -integral (|^) can be performed by the 
saddle point method. The maximumQ of 8 is reached for r = and ft, = 0. 
Expanding near this point to second order in r and 0j one finds 

© = 3^(AVM S )(1 - 2r 2 ) - ^C( Mg , Ms )^^|. (13) 

Finally, integration over r and 9i gives in the saddle point approximation the 
physical canonical partition function (after Wick rotation) 



ZQ(/3,v,fi q ,(x s ) = -V27rC (fi q ,fi s )u 2 ( — ) exp 



1 itw 

3^ 



(14) 



where u(fi q ,fi s ) are known El functions of quark chemical potentials and gluon 
and quark degrees of freedom (see Eqs. fl3l| ) and ( |32"| ) below for general case of 
N colors). 

3.2 The bag density of states and BBSBM mass spectrum 

The density of states <r(W, v, fi q , fx s ) of a bag of volume v, at energy W, is 
derived from its partition function Zq((3,v, fx q , /i s ), Eq. (|l4|), through inverse 



Laplace transform: 

a(W,v,n q ,fi s 



1 



C+200 



dfi exp{l3W)Z Q ((3 1 fi q , li s ) . 



(15) 



c+ioo 



For large W this integral can be evaluated by the steepest descent method to 
obtain the asymptotic bag level density: 

a(W, v, N , fx s ) = Ic-W^v^^W- 11 ' 2 exp U/3 u^v^W 3 / 4 ] . (16) 
8 L 

The GCPF of a gas of hadrons (bags) with mass m., and proper volume Vi 
can written in the Boltzmann approximation Q 



z H (T,v,\ g ,x s )=Y,w /n 

N=0 ' i=l 



N 



x e -W+^ T ( mi)Ui ,A g) A s ) 




(17) 



where T(m,v, X q , X s ) dmdv is the number of bag states in the mass interval 
(m, m + dm) and volume interval (w, w + dw). It has been shownB that the v de- 
pendence of r(m, w, A g , A s ) is obtained through the derivative of cr(W, v, X q , A s ), 
as given by Eq. (tlq) with respect to the volume v, 



d 

T(m,v,X, A s ) = — a(W,v, X, X s ) 
ov 



(18) 

W—m — Bv 

where the bag model relation W = m — Bv between the energy W and the 
mass m is used, B being the bag constant. This yields the BBSBM asymptotic 
hadronic mass spectrum: 

r(m, v, A, A s ) « V 5 / 4 i;- 9 / 4 (to - Bv)~ 19 ^ 



x exp 



3 v ; 



(19) 



This is the asymptotic expression of r(m, v, X q , X s ) derived from the asymptotic 
expression of cr(W, V, A, A s ). For small masses a phenomenological function To 
is introduced 



(m,v,Xg,X s ) = ^gi(X qi X s ) 8{m - m t ) 5(v - , 



(20) 



where gi(X qi X s ) is the multiplicity modified by fugacities of a hadron of mass 
to,- and volume Vi. 



4 Description of the HG — » QGP phase transition 

4-1 Pressure partition function 

We study here the analytic properties of the GCPF to show that the sys- 
tem undergoes a phase transition from hadron to QGP phase, and we imple- 
ment the generalization to N-colors. The analytic properties of the GCPF 
Zh(T, V, \ q , A s ), Eq. [I?], in the thermodynamic limit, (V — > oo), are most eas- 
ily investigated by considering the corresponding pressure-GCPF which is its 
Laplace transform with respect to the volume: 



/>oo 

U{T,s,X q ,X s ) = / dVe- sV Zv(T,V,\ q ,\ s ). 
Jo 

Substituting Eq. ( |l7| ) in Eq. (^l|) , one obtains 



(21) 



II(T, S ,A,A S ) = — - (22) 

s - f{P, s, X,X S ) 



where 



f(0,s,\ g ,\ s )=f o + 1 ^ 3 dv dmm 2 K 2 (m(3)T(m,v,---)e- sv , (23) 

'v J Mo+Bv 



2tt 2 /3 



and where /o = /o(/3, s, \ q , \ s ) corresponds to the contribution coming from 
To. The lower limits Mq and Vq have values above which the asymptotic 
expression (|l9|) begins to hold. 

Because of the properties of the inverse Laplace transform the behavior 
of Zh (T, V 1 X q , X s ) in the thermodynamic limit is governed_hy the rightmost 
singularity of IT(T, s, \ q , X s ). From Eq. ( |2^ ) it can be shownE3 that LT has two 
singularities on the real axis of the s-plane: a pole so given by 

8 (J3,\ 9 ,\ B ) = f(J3,a (J3,\ 9 ,\.)), (24) 

and a singularity s c (/3, A 9 ,A S ) induced by f(/3,s,X q ,X s ) which remains finite 
for s — > s c . Furthermore s c is found to be given by 

a c (/3, X q , X a ) = ^u(X q , X s )(p- 3 - BP) . (25) 

For given values of X q and A s the transition temperature Tc(X q , X s ) = is 
determined by matching of these two singularities, namely when 

s c ([3 c (X), X) = s (/3 c (A), A) = /(/3 C (A), s c (/3 c (A), A), A) , (26) 

where we have written A instead of (X q , X s ) to simplify the notation. 
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4-2 Properties of hadronic phases 



For T < Tc, in the hadron phase, the rightmost singularity is So(/3, A 9 , A s ) and 
the partition function is governed by this singularity in the thermodynamic 
limit. Then the pressure and the energy density are given by 



P H (T, A) = lim ^ In Z H (T, V, A) = Ts (fi,X), 
e H (T, A) = - lim — \nZ H (T, V, A) = — — s (P, A) . 

V^oo Op Op 



(27) 



For T > Tc(\ q , A s ), the leading singularity is s c and we have in the thermo- 
dynamic limit 



P Q (T,A)= Ts c (0,X) 



■ u(X) T 4 - B . 



e Q (T, A) = -QpS c (0, A) = u(X) T 4 + B . 



(28) 



Eqs. ( p8|) are the equations of a non interacting QGP. Of course other quan- 
tities like baryon number density, quark number densities can be derived form 
the expression of the partition function in the thermodynamic limit. Fig. [I] 
displays the energy density and the pressure for both phases for a value of 
B 1 ^ = 225 MeV, corresponding to T c = 150 MeV. Fig. | is the phase di- 
agram obtained by solving Eq. ( p6| ) for different values of the strange quark 
chemical potential fi s . 



4-3 Extension to N-colors 

The previous results can naturally be extended to SU C (N). By the same 
method as in Sec. 3.1 and Sec. 3.2, for N = 3, we obtain for the grand canonical 
partition function 



Z Q (P, v, N ^ s ) = A C ^- N2 ^u-^(^) N2+i)/2 exp 



1 uv 

3/?3 



(29) 



where 



N-l 



A = 



8(2tt 



,(JV-4)/2 



N 



(30) 



8 



\ 

u 



E- 




0.1 



0.15 



0.2 
T [GeV] 



I 

0.25 



0.3 



Figure 1: Energy density and pressure versus temperature: solid curve gives the EoS in the 
case fi q = [i a = and for B 1 / 4 = 225 MeV corresponding to T = 150 MeV. Dashed line is 
the pressure in the same conditions. 



u(jj, q ,n a ) = — 



15 N, fpn q 



(N 2 -l)d g + -N(d q + d s ) + —d q (^L) 1 



-—Nd s 
2 



7T 



2tt 2 



and 



AT d Q + ds d„ ( (3uo 
jd g + — — + - 



Put 



•2,, 2 

q 



P 2 V 



2vr 2 



(31) 



(32) 



By inverse Laplace transform of Zq([3, v, fi q , fj, s ), we obtain the bag density 
of levels a(W,v, /j, q , fi s ) for large W — m — Bv, from which the hadronic mass 
spectrum r(m,i>, fj, q , fj, s ) for large m is derived, according to Eq. ( |lg| ) 

r(m,«,M g ,M a )~C (1 -^ )/ V( 1 +^8 w 7^e^[4/3« 1 /V/ 4 W 3 / 4 ] , (33) 
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Figure 2: Critical curves in the plane T-fi q for fi s = 0, 100, 200, 300, 400, 500, 536, 560 
MeV, right to left. The lower curve corresponds to K condensation. 



where fcj=. -(N 2 + 9)/8 and 5 = -(11 + 3N 2 )/8. Furthermore it can be 
shown E2rllll that the conditions for a phase transition to occur are 7 + S = 
-(5 + N 2 )/2 < -3 and S = -(11 + 3iV 2 )/8 < -7/4. Both conditions give 
N 2 > 1. This can also be seen on the mass dependence of the spectrum. From 
the bag relations, rn = 4Bv = W + Bv, we have v = m/AB and W = 3m/ 4. 
Substituting for W and v in Eq. ( |33| ) we find 



T(TO ;/Ws ) ~ m -(^+ 5 )/2 c{ l~NZ)/2 u 3 { l + N*)/8 B (9 + NZ)/8 



exp 



( — ) 

\3BJ 



1/4 



m 



_ _ (34) 

where u and C are given by Eq. (pll) and Eq. ( |32| ) . Comparison with a spectrum 
of the form in Eq. (^) shows that a = —(N 2 + 5)/2 which is, for TV > 1, smaller 
than —7/2 as required for a phase transition to occur (see Sec. |2.2| ). Thus for 
N = 1 there is no phase transition. For N — 2 our approach is not-equivalent 
to lattice-QCD where one find a second order phase transition 113, a result 
requiring further study. 
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5 Conclusions 



Needless to repeat here, lattice QCD and its continuum limit is the fundamen- 
tal theory for describing in a unified scheme the deconfinement phase transi- 
tion, also for finite baryon density, and we are looking forward to future devel- 
opments in this direction. As we await these advances, we have offered here 
a relatively simple extension of the statistical bootstrap model, motivated by 
some well understood properties of strong interactions. We have shown that in 
theoretically consistent manner we can study the properties of deconfinement 
phase transition at finite baryon density. There are several simplifications in 
our approach which need attention: for example we assumed vanishing quark 
mass even for the s-quark, and we ignored QCD interactions in the bag. Even 
so, we find several interesting results: the properties of the latent heat of the 
phase transition, and the general features of the hadronic mass spectrum are 
in as one could wish for, and there is no phase transition when N = 1. 
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